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Utilizing a SU(2) gauge symmetry technique in the quasiclassical diffusive regime, we theoretically study 
finite-sized two-dimensional intrinsic spin-orbit coupled superconductor/normal-metal/superconductor (S/N/S) 
hybrid structures with a single spin-active interface. We consider intrinsic spin-orbit interactions (ISOls) that 
are confined within the N wire and absent in the s-wave superconducting electrodes (S). Using experimentally 
feasible parameters, we demonstrate that the coupling of the ISOIs and spin moment of the spin-active interface 
results in maximum singlet-triplet conversion and accumulation of spin current density at the comers of the N 
wire nearest the spin-active interface. By solely modulating the superconducting phase difference, we show 
how the opposing parities of the charge and spin currents provide an effective venue to experimentally examine 
pure edge spin currents not accompanied by charge currents. These effects occur in the absence of externally 
imposed fields, and moreover are insensitive to the arbitrary orientations of the interface spin moment. The 
experimental implementation of these robust edge phenomena are also discussed. 

PACS numbers: 74.50.+r, 74.25.Ha, 74.78.Na, 74.50.+r, 74.45.+C 


I. INTRODUCTION 

The interaction of a moving particles’ spin with its lin¬ 
ear momentum embodies the so-called spin-orbit interaction 
(SOI). The SOI is a quantum mechanical effect that is rela¬ 
tivistic in origin. For materials possessing a strong SOI effect, 
it becomes possible to manipulate spin currents with less dis¬ 
sipation, higher speeds, and lower power consumption com¬ 
pared to conventional charge-based devices*. Consequently, 
a number of high-performance devices that exploit the SOI 
effect have been proposed, including, spin transistors, and 
devices that store or transport information^"^. The types of 
SOIs can be categorized as follows: i) intrinsic (originating 
from the electronic band structure of the material) and ii) ex¬ 
trinsic (originating from the spin-dependent scattering from 
impurities)^. Of particular interest are intrinsic SOIs (ISOIs), 
due to their controllability by tuning a gate voltage"*’®’^"*^. 
Two commonly studied ISOIs are the Rashba and Dresselhaus 
types. The Rashba SOI'^ can be described via spatial inversion 
asymmetries, while the Dresselhaus SOI*^ is a result of bulk 
inversion asymmetries within the crystal stnicture"*’®. 

There have also been extensive efforts to manipulate 
the spin currents*’^’^’*"*’*^ in SOI systems via the spin Hall 
effect*®"^^, and the quantum spin Hall effect^^’^^. Since spin 
curi'ents are weakly sensitive to nonmagnetic impurities and 
temperature^, more opportunities arise in the development of 
high speed low-dissipative spintronic devices*. Along these 
lines, superconducting heterostructures have been making 
strides as potential platforms where spin-orbit coupling (SOC) 
plays a key role, including scenarios involving the spin-Hall 
effect When considering superconducting hybrids 

with SOC, interface phenomena at superconducting junctions 
becomes particularly important. For example, the interface 
of a hybrid superconducting junction can behave as a spin- 
polarizer when it is coated by an ultrathin uniform F layer. 
The study of interface effects that involve spin-dependent 


scattering^**"^^ has spanned considerable theoretical'***’"**’^^"^* 
and experimental works'*®"^*. Advancements in nanofabrica- 
tion and theoretical techniques involving superconducting hy¬ 
brids with spin-active interfaces have thus created new venues 
for controlling superconducting pair correlations, spin cur¬ 
rents, and majorana fermions"***^^’^^. 

To explore the interplay of these phenomena, we consider 
charge and spin currents in a finite sized intrinsic spin-orbit 
coupled s-wave superconductor/normal-metal/s-wave super¬ 
conductor (SINIS) junction with a single spin-active in¬ 
terface. We utilize a spin-parameterized two-dimensional 
Keldysh-Usadel technique^® in the presence of ISOIs. In or¬ 
der to theoretically account for spin-polarization and spin- 
dependent phase shifts that a quasiparticle experiences upon 
transmitting across spin-active interfaces, spin-boundary con¬ 
ditions are utilized^^. The spin-parametrization scheme al¬ 
lows us to isolate the spin-singlet and spin-triplet cori'elations, 
and pinpoint their spatial behavior. We find that the combi¬ 
nation of interface spin moment and ISOIs results in triplet 
pairings with m = 0, ±1 spin projections along the quantiza¬ 
tion axis^^’^®’^^"^^. We also find that maximum singlet-triplet 
conversion and spin-current densities takes place at the cor¬ 
ners of the N wire nearest the spin-active interface, where 
the spin accumulation is greatest. The spin currents possess 
three nonzero spin components, independent of either the ac¬ 
tual type of ISOI present in the N wire or spin moment orien¬ 
tation of the S/N spin-active interface. When comparing the 
spin and charge currents as functions of the superconducting 
phase difference, ip, we show that current phase relations for 
the charge supercurrents are typically governed by sinusoidal- 
like, odd functions in (p, although anomalous behavior^**’^*’®* 
can arise. The spin currents however, are even functions of 
These opposing behaviors of the spin and charge cur¬ 
rents present a simple and experimentally feasible platform 
to effectively generate proximity-induced spin-triplet super¬ 
conducting pairings and edge spin currents in the absence of 
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any charge supercurrent. It was demonstrated in Ref. 59 that 
the combination of spontaneously broken time-reversal sym¬ 
metry and lack of inversion symmetry can result in sponta¬ 
neously accumulated spin currents at the edges of finite-size 
two-dimensional magnetic S/F hybrids. Moreover, we de¬ 
scribe experimentally accessible signatures in the physically 
relevant quantities, and discuss realistic material parameters 
and geometrical configurations that lead to the edge spin ef¬ 
fects predicted here. Our proposed hybrid structure, based 
on its intrinsic properties alone, can be viewed as a sim¬ 
pler alternative to differing systems that rely inextricably on 
externally imposed fields to generate the desired edge spin 

cuiTentsi®'23,26_ 

The paper is organized as follows: We outline the theoret¬ 
ical techniques and approximations used to characterize the 
intrinsic spin-orbit coupled superconducting S/N/S hybrid 
structures with spin-active interfaces in Sec. II. In Sec. Ill, 
we discuss all of the types of superconducting pairings present 
(spin-singlet and spin-triplet correlations) and illustrate the as¬ 
sociated spatial profiles, which follow directly from the inher¬ 
ent proximity effects. Next, we present results for the spin 
current densities, with spatial maps, and discuss possible ex¬ 
perimental realizations of the presented edge spin phenomena. 
In addition, we discuss the spin current symmetries relative to 
the spin moment orientation of the spin-active interface in the 
presence of Rashba or Dresselhaus SOC. We finally summa¬ 
rize our findings in Sec. IV. 


II. THEORETICAL FORMALISM 

A quasiclassical framework has recently been developed for 
superconducting hybrid structures in the presence of generic 
spin-dependent fields^®. These generic fields can be reduced 
to ISOIs, such as, Rashba®^ and Dresselhaus^^, in terms of the 
quasiparticles’ linear momentum [p = {px,Py,Pz)]- If we 
define a vector of Pauli matrices r = {Tx,Ty,Tz) (see Ap¬ 
pendix A), the corresponding Hamiltonians describing these 
ISOIs can be straightforwardly expressed as, 

"Hr = Or(p Xt) ■ Z ^RiPxT"^ - PyT""), 

■Hd = ^d{P ■ t) OD(pa:T“ -f PyT^), 

where the momentum is restricted to the xy plane. The R 
and D indices represent the Rashba and Dresselhaus SOIs 
with strength Hr and Ho, respectively. A linearized ISOI 
can be treated as an effective background field that obeys the 
SU(2) gauge symmetries.Therefore, to incorporate 
ISOIs into the quasiclassical approach, it is sufficient for par¬ 
tial derivatives to be interchanged with their covariants. 

This simple prescription is one of the advantages of the SU(2) 
approach, besides the convenient definitions of physical quan¬ 
tities such as the spin currents. 

A description of quasiparticle transport inside a supercon¬ 
ducting medium is provided by the Dyson equation.^* The 
corresponding equation of motion in the quasiclassical ap¬ 
proximation for clean systems reduces to the so-called Eilen- 
berger equation^'. The Eilenberger equation can be further 
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FIG. 1. (Color online) Schematic of a two-dimensional SIN IS 
Josephson junction with a spin-active interface at x = 0 described by 
S; = (S'f, S'f, Sf ). The intrinsic spin-orbit coupled normal-metal 
strip (A) is of length and width djv, and Wjv, respectively. The su¬ 
perconducting electrodes (Sc), however, are SOTfree and connected 
to the N wire at x = 0, and x = d^- The two-dimensional junction 
resides in the xy plane so that the SIN interfaces are parallel with 
the y axis. The cone depicts the tip of a scanning tunneling micro¬ 
scope (STM) that can sweep the entire surface of the N layer in the 
xy plane. 


reduced to a simpler set of equations in the diffusive regime, 
where the quasiparticles are scattered into random directions. 
This permits integration of the Eilenberger equation over all 
possible momentum directions, yielding a simpler picture for 
highly impure systems, as first introduced by Usadel^^. The 
resultant Usadel equation is the central equation used in this 
paper, and can be expressed compactly 




d,G{r,e)[d,G{r,e)\ +i[ep 3 


A,G(r,e)] = 0, (1) 


where pa is a 4 x 4 Pauli matrix (see Appendix A), S) repre¬ 
sents the diffusion constant, and A is a 4 x 4 matrix which 
represents the superconducting gap^^. We denote the quasi¬ 
particle energy by e, relative to the Eermi energy, ep- In the 
superconducting leads, the Usadel equation, Eq. (1), is solved 
in the presence of A which results in the BCS bulk solution 
given by Eq. (6). Within the nonsuperconducting region how¬ 
ever, A = 0, and the boundary conditions, Eq. (4), should 
be simultaneously satisfied. The total Green’s function for the 
system, G(r,e), is comprised of the Advanced (G^(r,e)), 
Retarded (G^(r, e)), and Keldysh (G^ {r, e)) propagators: 


G(r,e) = 
G^{r,e) = 


G^{r,e) G^{r,e) 

0 G^(r,e)J’ 

-G{r,-e) -F{r,-e) 
F*ir,e) g*ir,e) 


Since we consider the low proximity limit of the diffusive 
regime,^^, the normal and anomalous components of the 
Green’s function can be approximated by, Q{r,e) ~ 1 and 
F{r,£) <C 1, respectively. Thus the advanced component of 
total Green’s function reduces to: 


G^{r,e) 


-1 -F{r,-e) 

F*{r,s) 1 


( 2 ) 


Within the low proximity approximation, the advanced com- 
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ponent can be ultimately expressed as; 


conducting bulk solution, G§q^, by piecewise functions: 


G^ir,e) = 

-1 0 -fn{r,-e) 

0 -1 -/+(r,-£) 

f^{r,s) flir,e) 1 

_/+(’’,£) 0 

In equilibrium, the Retarded and Keldysh blocks 
of the total Green’s function are obtained via: 
G^{r,e) = -{p 3 G^{r,e)p 3 y, and G^{r,e) = 

tanh{ekBT/2){G^{r,e) — G^{r,e)'\. The Boltzmann 
constant and system temperature are denoted hy ks and T, 
respectively. Within the low proximity approximation, the 
linearized Usadel equation involves sixteen coupled complex 
partial differential equations which become even more 
complicated in the presence of ISOI terms. Unfortunately, 
the resultant system of coupled differential equations can be 
simplified and decoupled only in extreme limits that can be 
experimentally unrealistic.When such simplihcations 
are made, the equations lead to analytical expressions for the 
Green’s function components'^. For the complicated system 
considered in this paper however, computational methods 
are the most efficient, and sometimes only possible routes to 
investigate experimentally accessible transport properties.^® 
The complex partial differential equations must be supple¬ 
mented by the appropriate boundary conditions to properly 
describe the spin and charge currents in spin-orbit coupled 
SIN IS hybrids with spin-active interfaces;^^’®'^ 


0 

1 


■ ( 3 ) 


s(£) = sinh ■!?(£) = 


-A 


sgn(£) 


■\J 

c{e) = coshi?(£) = 

I £' 


0(£" - A") - 




,0(A^-£^) , 


'J 


0(£" - A") - 


IS 




= 0(A^-£^), 


where 0(x) is the Heaviside step function and A is the super¬ 
conducting gap at temperature T for a conventional s-wave 
superconductor. 

We next employ a spin-dependent field technique that 
permits the incorporation of ISOIs into the Keldysh-Usadel 
approach^®. As stated earlier, this technique has been widely 
used in the literature^and was most recently ex¬ 
tended for superconducting heterostructures^®. In much the 
same spirit, we adopt a generic tensor vector potential A{r) — 
{A^ {r),Ay{r),A^ .39,59,67,69 


Ax{r) = 


''+Ayir)ry+AUr)T^}, 

(7a) 

^y{r) = 


^'+Ayy{r)Ty+A;{r)T^}, 

(7b) 

A^{r) = 


'+A^(r)r^+A^(r)T"}. 

(7c) 

We can now dehne the covariant derivative by. 



d = 

dl — ieA{r), 

(8) 


C,{G{r,£)dG{r,£)) ■ n = [GBcs{Oir,£),G{r,£)] 

±i[Sir ■ v>,G{r,£)], (4) 


where the unit vector, n, is directed normal to an interface, 
and it is assumed for the time being that the left and right 
interfaces [Fig. 1] are spin-active. The leakage intensity of 
superconducting correlations from the S electrodes to the N 
wire is controlled by the ratio between the resistance of the 
barrier region Rb and the resistance in the normal layer 

c = Rb/Rn.^' 

We describe the spin moments of the left (1) and right (r) 
interfaces by two generic vectors as follows:"^®’'*' 

Sir = {SI. SI, SI), ( 5 ) 

where the Sir can have arbitrary directions and magnitude of 
the spin moment at the SIN interfaces. The solution to Eq. (1) 
for a bulk even-frequency s-wave superconductor results in 


G'bCs(^i 


(T® cosh "&{£) sinh '!?(£) 

icrJ^e-*®sinh^l(£) -ct® cosh^l(£) J ’ 


( 6 ) 


where d = {dx,dy,dz). Hence, the brackets in the Usadel 
equation Eq. (1) and the boundary conditions Eq. (4) (as well 
as the charge and spin currents discussed below) are equiva¬ 
lent to: 


[r?, G(r)] = dG{r) — ie[A{r),G{r)\. (9) 

It should be noted that the quasiclassical approach employed 
here allows for the study of systems with spin dependent vec¬ 
tor potentials possessing arbitrary spatial patterns, and spin- 
active interfaces with arbitrary spin moment directions. Here 
we assume that the impurity scattering (encapsulated by the 
diffusion constant) is spin-independent, and thus the spin- 
dependent helds introduced in Eqs. (7) describe the spin-orbit 
coupling for the system. Within the quasiclassical approxi¬ 
mation, the quasiparticles’ momentum is localized around the 
Eermi level. Therefore, the spin moment amplitude for spin- 
active interfaces \Spr\, the vector potential, |A|, and super¬ 
conducting gap |A|, should all be appropriately smaller than 
the Eermi energy £j’®®. 

A specihc choice of the tensor vector potential®^ [Eq. (7)] 
that results in linearized Rashba (a ^ 0, /3 = 0)®^ and Dres- 
selhaus (/3 ^ 0, a = 0)'^ SOCs is: 


where !){£) = arctanh(|A|/£), and 6 represents the macro¬ 
scopic phase of the bulk superconductor. The phase differ¬ 
ence between the left and right S electrodes, shown in Eig. 1, 
is denoted 6i — 9r = <p- We dehne s and c terms in the super- 


= -Al = 2/3, 

Al = -Al = 2a, 

Al = Al = 0, 

— AX _ Ay _ n 
■^z *^ 2 : ^5 


( 10 ) 
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where a and /3 are constants and determine the strength of 
Rashba and Dresselhaus SOIs. This choice simplifies the re¬ 
sultant Usadel equations since da = 9/3 = 0^ d-A{r) = 0. 
Hence, by substituting the above set of parameters, Eqs. (10), 
into Eqs. (7) we arrive at:^® 


=/3t“ - ar^, (11a) 

Ay = ar^ — . (11b) 


Although these assumptions lead to further simplifications of 
the Usadel equations, the end result involves sixteen coupled 
complex partial differential equations that we have analyti¬ 
cally derived, but omitted here due to their excessive size. 

As mentioned in the introduction, crystallographic in¬ 
version asymmetries^^ or lack of structural inversion 
symmetries®’®®’^^’^* in heterostructures may cause the finite 
ISOIs considered in this paper. Eor example, engineered strain 
can induce such inversion asymmetries^®’®^’®^’^^ and conse¬ 
quently, ISOIs. Alternatively, the adjoining of differing ma¬ 
terials may generate interfacial Nonetheless, 

currently there is no straightforward method to measure SOIs 
in a hybrid structure. One approach might involve first princi¬ 
ple calculations for combined materials. Also, photoemission 
spectroscopy*® and spin transfer torque experiments can pro¬ 
vide realistic values for the ISOIs. 

One of the most striking topics in the study of transport 
in junction systems involves spin currents. Since a decade 
ago, various features and behaviors of spin currents in hybrid 
structures have extensively been studied.The 
spin and charge currents are key quantities that reveal useful 
insights into the system transport characteristics. These phys¬ 
ical quantities are also crucial to the application of nanoscale 
elements in superconducting spintronics devices.*®’^* The vec¬ 
tor charge iJA and spin (J^'*') current densities can be cal¬ 
culated by the Keldysh block {K) when the system is in 
equilibrium;®^’®® 




^+oo 

deTr|p3(G(r)[a,G(r)])^| 


, ( 12 ) 




r+oo 

deTx 


{p3i"®'(G(r)[9,G(r)])^|| 


(13) 


where Jq = NqbD/A, Jq = hjQ/2e, and Nq is the number 
of states at the Eermi surface. The vector current densities 
provide local directions and amplitudes for the currents as a 
function of position. We designate 7 = x,y,z for the three 
components of spin current, e.g., represents the x com¬ 
ponent of spin current. The integral of Eq. (12) over the y 
direction, shown in Eig. 1, provides the total charge supercur¬ 
rent flowing across the system. 

In order to pinpoint the behavior of the spin-singlet and 
spin-triplet pairings inside the spin-orbit coupled N wire, 
we exploit a spin-parametrization scheme*®’®*’®"^ where the 
anomalous component of the Green’s function, Eq. (2), can 
be parameterized as follows: 


=S(r,£)+cr-T(r,e) = 

cr°S(r, e) + a'-‘Ta;{r, e) + (J^Ty{r, e) + cr^T^jr, £). (14) 


In terms of this spin-parametrization, the quantities S and T^^y 
then correspond to the singlet and triplet components with 
m = ±1 total spin projection along the spin quantization axis, 
while Tz represents the triplet component with m = 0 .*®’^'’®* 
Here, the spin quantization axis is assumed fixed along the 
z direction throughout the entire system. In a uniform ferro¬ 
magnetic system, it has been demonstrated that the singlet and 
m = 0 triplet component decay rapidly while the m = ±1 
triplet components, if exist, can propagate over longer dis¬ 
tances compared to the former correlations. By noting this 
aspect of triplet correlations, namely, the degree of their pen¬ 
etration (i.e. the distance that the correlations are nonzero) 
into a system, one may classify them as ‘short-range’ and 
‘long-range’ correlations. Considering this classification, the 
m = 0 triplet component in a uniform ferromagnet is short- 
ranged while the m = ±1 are long-ranged.®® ®® Therefore, the 
parametrization scheme we utilize allows for explicit determi¬ 
nation of the spatial profiles for the different superconducting 
pairings in intrinsically spin-orbit S/N/S systems with spin- 
active interfaces and hence, their short-range and long-range 
natures. 


III. RESULTS AND DISCUSSIONS 


For additional insight and comparison purposes, we first 
consider a S/N/S junction with negligible SOIs and either 
one or two spin-active interfaces. We compute the spin cur¬ 
rents, and discuss the singlet and triplet correlations in such 
systems. Next, we compute these same quantities after incor¬ 
porating one of the Rashba or Dresselhaus ISOIs introduced 
above. As remarked earlier, Eq. (1) in the presence of ISOI 
terms leads to lengthy coupled partial differential equations. 
Although we have analytically derived the current densities 
[Eqs. (12) and (13)] for numerical implementation, they lead 
to cumbersome expressions that are not repeated here. In what 
follows, we normalize the quasiparticle energy £ by the su¬ 
perconducting gap at zero temperature Aq = A(T = 0 ), 
and all lengths by the superconducting coherence length 
which is defined as yJfiT)/N q. The barrier resistance at the 
S/N interfaces is set to = 4. This value of the barrier re¬ 
sistance warrants the validity of the low proximity limit i.e. 
Q{r,e) ~ 1 and /F{r,e) <C 1. We consider a fixed value 
for the spin moment amplitude of the spin-active interfaces, 
= GT{4>ir)~^ = 15, corresponding to realistic experi¬ 
mental situations.®® Here, Gt is the barrier conductance and 
(pir represents the spin-dependent interfacial phase-shifts at 
the left and right interfaces .^®’‘^®’‘^® We note that by considering 
other values for | Sir\, within reasonable experimental bounds, 
the influence on the results is negligible. In our actual com¬ 
putations, we have adopted natural units, so that = 1 . 

To find stable solutions to Eq. (1), and thus to obtain currents 
given by Eqs. (12) and (13), we have added a small imaginary 
part, S ~ O.OIAq, to the quasiparticles’ energy e ^ e + i6. 
Physically, the additional imaginary part can be considered as 
a contribution from inelastic scatterings.®* Due to this imagi¬ 
nary part, we have taken the modulus of quantities (denoted by 
the usual |...|). Also, for symmetry reasons we restrict the spa- 
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FIG. 2. (Color online) Spatial profiles of the three spin current components, y), y), and J^^{x, y) in a S/N/S junction with 

spin-active interfaces and no spin orbit interactions. Top row panels are against x (positions along the junction length) at differing locations 
along y = 0.255s, 0.55s, 0.755s, l-05s (see Fig. 1). By contrast, the bottom row panels are functions of y, plotted at four positions along 
the junction length x = 0.255s, 0.55s, 0.755s, l-05s- For the N wire we have, (In = Wn = 2.05s, and the phase difference between the S 
electrodes is fixed at = k/2. The spin moments of the spin-active interfaces are denoted by Sw = {Sf^, Sf^., Si^) at the left(Z) and right(r) 
interfaces, (a) The right interface is spin-inactive while the spin moment of left interface points along the z direction. Si = (0,0, iS'f ). In 
(b), the two interfaces are spin-active with parallel spin moments Sir = (0, 0, while in (c). Si = (0, 0, iSf ) and Sr = (0, ±5;!, 0). 


tial profiles to the regions 0 < a: < dj^/^ and 0 < y < Wn/2 
in the figures presented throughout the paper. 


III.l. S/N/S junctions with spin active interfaces in the 
absence of SOC 

Figure 2 displays the spatial profile of the spin current com¬ 
ponents in a S/N/S junction containing spin-active interfaces 
and negligible SOIs. As seen in Fig. 1, the junction resides 
in the xy plane so that the S/N interfaces are located at 
X = 0 , djv and the vacuum borders in the y direction are found 
at y = 0, Wn- The junction width Wn and length dAr are set 
to a representative value of 2 . 05 s, without loss of generality. 
The top set of panels in Fig. 2 are functions of the x coor¬ 
dinate at differing y, i.e., 0.255s, 0.55s, 0.755s, and 1.05s, 

while the bottom row of panels are functions of y at fixed x = 

0.255s, 0.55s, 0.755s, and 1.05s- In Fig- 2 (a), the spin mo¬ 
ment of the left interface is fixed along z. Si = (0, 0, iSf), 
while the right interface is spin-inactive, Sr = 0. Note that 
we have denoted vector S by three ‘positive’ scalar entries 
S^, S'^^ S^, representing projection of the vector in the x, y, z 
directions (in the Cartesian coordinate). Here, thus, ± signs 
indicate the orientation of that component which can be par¬ 
allel (H-) or antiparallel (-) to specific directions x, y, z. It can 
be seen that the only nonzero component of spin current is 
J^^{x,y) due to Si = (0,0,±S'f). This component of spin 
current is maximum at the left S/N interface, x = 0, because 
of the abrupt spin-imbalance, and drops to a vanishingly small 
value near the right interface, x = dN, which decays to zero 
over the much longer dN length scale. Moreover, the plot 


demonstrates a uniform distribution of spin current along y. 
We see that the curves at various y locations overlap, consis¬ 
tent with the bottom panel of Fig. 2(a), where the spin current 
is constant for a given x. 

As described above and seen in Fig. 2, the results for the 
two-dimensional S/N/S junction in the absence of ISOIs re¬ 
duces it to an effectively one-dimensional problem. Hence, 
to gain more insight, we derive analytical expressions for the 
charge and spin current densities in a simple structure, namely 
a one-dimensional S/N/S junction with a single spin-active 
interface. To this end, we first derive solutions to the com¬ 
ponents of the Green’s function, Eq. (3), i.e. f±{r,zLe), 
[f±{r,±e)]*, and ie)]* in a one¬ 

dimensional S/N/S junction where Si = (0,0,±S'f) and 
Sr = 0. If we define a normalized coordinate x = x/dN, 
we end up with the following expression for /-(x, —e) (sim¬ 
ilar expressions arise for the other components): 


/-(x,-e) = s*(e) X 

-iip /2 C^'lcosh Ax + cosh A(1 — x)] — iSf sinh Ax 
5-^ [C-^ sinh X — iSf cosh A] 


To simplify the solutions, we denote 2iejeT = A (in which 
St is the Thouless energy) and assume Oir = Supj^. By sub¬ 
stituting the obtained solutions into Eq. (12) we arrive at the 
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following expression for the charge current: 

Jxix,^) = JoJ d£2iXtanh{^^^)'^smip, (15) 
= 

[s*(-e)]^ sin Ajs'f^ - + (S'f^ + C^A^) cosh2A} 

+[s*(e)]%inhA|5'fVC^A2 + - C^A^) cos2a}, 

= 

|s'f^ cos^ A + C^A^ sin^ cosh^ A + C^A^ sinh^ a|. 

It is immediately evident that the charge current has the usual 
sin ip odd-functionality in the superconducting phase differ¬ 
ence. The same procedure is followed to derive analytical 
expressions for the spin current components using Eq. (13). 
Since S'; = (0,0, ±Sf), we find </?) = J^y{x,p) = 0 

and (fi) ^ 0. Unfortunately, even by means of the sim¬ 

plifying approximations made to the equations thus far, we 
arrive at a rather lengthy expression for the z component of 
the spin current, J^^{x, Nevertheless, if we restrict our 
attention to the edge of the wire, x = 0, this spin current com¬ 
ponent reduces to the following: 




Mr 


Mr 


vr 


V 


sz 

2 






m 




^2 


},( 16 ) 


^ V 2 

[s*(—e)]^| 2 Sf^[l + cosAcos p] + C^A^[cos4A — 
6 cos A cos - 3] -f 2[Sf ^ -f -f {Sf + 
cos A cos p] cos2a|. 


[s*(£)]^|3[Sf^ -I- C^A^] coshAcos:/? -I- 
4Sf ^ cosh^ A + (Sf ^ — 3C^A^) cosh 3A cos p 
— 4C^A^[2 -I- cosh 2A] sinh^ a|, 

{Sf^ -f C^A^ -f [Sf - C^A^] cos2a}^ 

{Sf2 - C^A^ -f [Sf -f C^A^] cosh2A}^ 


The z component of spin current, is evidently an even- 
function of p, namely cos p, although it involves some com¬ 
plicated prefactors.This finding is consistent with Ref. 61. 
Note that when the other components of spin current, are 
present, the even-functionality in p holds, even in the presence 
of ISOIs. This effect is discussed further below. 

In Fig. 2(b), the right interface is now spin-active (at x = 
djsi). The spin moment direction of the spin-active interface at 
a: = 0 is intact while Sr = (0,0, ±5'^). As seen in the bottom 
panel, J‘^^{x,y) is still the only nonzero spin current com¬ 
ponent, which is constant in the y direction. The right spin- 
active interface causes an increase in J^^{x,y) at x = d^ 
due to a spin-imbalance effect. In Fig. 2(c), the spin moment 
of the interface at x = dAr is oriented along the y direction, 
i.e., Sr — (0, ±5'*(,0). We see that J®^(x,j/) and J^^{x,y) 
are both nonzero since Si and Sr are orthogonal. The spin 
current vanishes at the middle of the junction (x = dM/2.) 


and apparently the behavior of the components become in¬ 
terchanged at this location. From the bottom row of panels 
in Fig. 2, we see that the spin-active interfaces with various 
spin moment orientations would lead to uniformly distributed 
spin currents along the junction width in the y direction. In 
other words, the spin-active interfaces alone are unable to 
induce any spin accumulation at the vacuum borders of the 
N vvire.'^’'®”^^’®^’® The triplet correlations in superconduct¬ 
ing hybrids with spin-active interfaces have extensively been 
studied."^®’'*^’"^^ In S/N/S systems with a single spin-active in¬ 
terface, no equal-spin pairing can arise (since a single quan¬ 
tization axis exists throughout the whole system), although 
opposite-spin triplets, T^(x,t/), can be induced. Figure 2(a), 
where only Jr is nonvanishing, confirms this phenomenon. 


III.2. Intrinsic spin orbit coupled S/N/S junctions with a 

single spin active interface: Singlet and triplet correlations 

Now we incorporate ISOIs in the S/N/S junction with one 
spin-active interface at x = 0 (see Fig. 1). Figure 3 ex¬ 
hibits first spatial profiles of the singlet (S(e, x, y)) and triplet 
(Ta;,y, 2 {e, X, j/)) correlations. Here we set S'; = (0, 0, ±Sf), 
Sr = 0, and assume that the ISOI is of the Rashba form, i.e., 
a ^ 0, and /3 = 0 (we later discuss the results of a Dres- 
selhaus SOC). We choose a = 2.0^^ as a representative 
value and emphasize that this specific choice has no influ¬ 
ence on the generality of our findings. In the low proximity 
limit, quasiparticles with low energies (e <C 1) tend to have 
the main contributions to the pair correlations. Accordingly, 
we therefore choose a representative value for the quasiparti¬ 
cles’ energy equal to e = O.lAg. The other parameters are 
kept unchanged. As clearly seen, the combination of an ISOI 
and spin moment of only one spin-active interface results in 
three nonzero components of the triplet correlations, which is 
in contrast to the case with zero ISOI shown in Fig. 2. This 
phenomena directly follows from the fact that the quasiparti¬ 
cle spin is tied to its momentum in the presence of an ISOI.^® 
The singlet |S(£,x,t/o)| is minimum at x = 0 and increases 
monotonically towards x = c?Ar at every point along the junc¬ 
tion width yo = 0.25^5,0.5^5,0.75^5, and 1 . 0 ^ 5 . This be¬ 
havior can be understood by noting the combination of in¬ 
terface spin moment and ISOIs at x = 0 abruptly converts 
the singlet correlations into triplet correlations. This picture 
is however reversed at x = djv where an ingredient to the 
singlet-triplet conversion is lacking, i.e., = 0 at x = dN- 

Examining the spatial map of the triplet correlations in Fig. 3, 
we find that the triplet correlations behave oppositely to the 
singlets. The triplet correlations T^, ^ ^(e, x, j/o) maxi¬ 
mum near x = 0, where the S/N interface is spin-active. 
Here, as remarked above, the combination of interface spin 
moment and ISOIs effectively converts the singlet supercon¬ 
ducting correlations into the triplet ones at x = 0. The triplet 
correlations decline as a function of x, and eventually convert 
into the singlets at x = c?Ar (at the spin-inactive interface). 
The triplet correlations, Tx,y, have nonzero spin-projections 
along the spin-quantization axis (m = ± 1 ) while m = 0 for 
the T 2 component. It is evident that T^ is drastically sup- 







7 




FIG. 3. (Color online) The spatial behaviors of the singlet S(e, x, y) and triplet T(£, x, y) correlations inside a Rashba S/N/S system with 
S; = (0,0, iS'f) and Sr = 0. The quasiparticles’ energy is set at e = O.lAo, djv = Wn = 2.0^s, and ip = tx 12. The top panels display the 
various pair correlations vs a; at y = 0.25^s, O.S^s, 0.75^s, and l.O^s, whereas the bottom panels exhibit the same quantities as a function of 
y at a: = 0.25^3, O.S^s, 0.75^3, and 1.0^3. 


pressed when moving away from the spin-active SIN inter¬ 
face at a; = 0 (a consequence of the so-called short-range 
behavior of T^). The spin-1 quantities, T^^y, on the other 
hand, remain nonzero over greater distances (the so-called 
long-range behavior). The actual distances that the triplet cor¬ 
relations can propagate over them before fully vanishing in 
a system depends on the system parameters such as temper¬ 
ature, degree of the interface opacity strength of the in¬ 
terface spin-activity and the magnitude of SOCs present in 
the system. Nonetheless, a direct comparison of penetration 
depth between the triplet correlations with zero total spin T^ 
and nonzero total spin i.e. llx,y clearly reveals that T^ is the 
short-ranged triplet component while Tlx,y are long-ranged in 
the ISO coupled SINIS junction with one spin-active inter¬ 
face. We note that this conclusion generally holds, indepen¬ 
dent of the representative values chosen. The bottom set of 
panels in Fig. 3 illustrates the superconducting correlations 
|S{£, xq, 2 /)|, and |Ta;_y_^(£, a;o, y)| as functions of y-position 
along the junction width where Xq = 0.25^s, 0.5^5,0.75^S> 
and l.O^s. The spatial distribution of the singlet correlations 
along y are unaffected by the coupling of the ISOIs and in¬ 
terface spin moment at a: = 0. The singlets, |S(£, xq, y)\, are 
constant along y, implying a uniform distribution along the 
junction width. This however differs considerably from the 
spatial behavior of the triplet correlations: The three triplet 
components x,y,z{£, xq, y)\ demonstrate appreciable accu¬ 
mulation at the transverse vacuum boundaries of the N wire 
(at y = 0, and y = Wn). Also, the results reveal that the max¬ 
imum singlet-triplet conversions occur at the corners of the N 
strip near the spin-active interface (near a; = 0, y = 0 and 
X = 0, y = Wn)-^'^ Further investigations have demonstrated 
that the maximum singlet-triplet conversion in such systems 
generally takes place at the corners of the N wire near any 
spin-active interfaces.^® We note that this finding is generic, 
robust, and independent of either interface spin moment di¬ 
rection or the actual type of ISOI considered. Similar spa¬ 
tial profiles appear when a = 0, and /3 0, or equivalently 


when a Dresselhaus SOI is considered. Our numerical inves¬ 
tigations have found that the corresponding Dresselhaus SOI 
results can be straightforwardly obtained by making the fol¬ 
lowing replacements in Fig. 3: |S’^(£, a:, y)| = |S°(£, a;, y)|, 
|Tf(£,a:,y)| = |T°(£,a:,y)|, |T^(£,a;,y)| = |T°(£,a;,y)], 
and |T^(£,a:,y)| = |T°(£, a:, y)|. The symmetries can be 
easily understood by considering the symmetries of the spin- 
dependent fields in Eq. (11), resulting in Rashba and Dressel¬ 
haus SOIs.^® The spin current components also show similar 
symmetries and we shall discuss them in detail at the end of 
this section. 


III.3. Intrinsic spin orbit coupled S/N/S junctions with a 
single spin active interface: Spin currents 

Figure 4 exhibits the corresponding spatial profiles of the 
spin current components, given by Eq. (13) for a Rashba 
S/N/S junction with one spin-active interface. We assume 
that the left interface at x = 0 is spin-active (see Eig. 1), 
and its spin moment is oriented along the y direction, namely. 
Si = (0, iS'f, 0) (and thus Sr = 0). To be specific, we first 
present the results of a Rashba (a ^ 0, /3 = 0) S/N/S system 
in our plots and then later expand our discussion to differing 
orientations of Si in the presence of either Rashba or Dres¬ 
selhaus SOIs. The top row of panels in Eig. 4 display the 
three components of spin current density Jp'{x,yo) flowing 
along X at yo = 0.25^s, O.S^g, 0.75^s, and l.O^s. The spatial 
variations of the x component provides the clearest and most 
useful information on the spin current behavior in such sys¬ 
tems. Therefore, we only present the x components in Eig. 4, 
while later discussing the vector spin current densities when 
presenting symmetries among the spin current components. 
The spin current density components vanish within x > cIn/^ 
while they are largest near the spin-active interface at x = 0. 
This behavior is consistent with the associated triplet correla- 
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S=(0,SJ',0), S =(0,0,0), a^O, ^5=0 (RashbaSOI) 



FIG. 4. (Color online) The spin current components, y), y), J^^{x, y) in a Rashba (a 7^ 0, /3 = 0) S/N/S junction with one 

spin-active interface (at x = 0). The interface spin moment is oriented along the y direction Si = (0, iS'f, 0) and Sr = 0 (see Fig. 1). The 
top panels show the spin currents vs a; at y = 0.25^s, O.S^s, O.FS^s, and l.O^s. The bottom panels exhibit the same quantities as a function 
of y at a; = 0.25^5, 0.5^s, 0.75^s, and l.O^s. The superconducting phase difference is set at y? = 7r/2, and (In = Wn = 2.0^s. 


tions investigated in the top row of panels in Fig. 3. The spin 
currents are zero dlx = cLm, where the SIN interface is spin- 
inactive. This hnding is also consistent with previous theoreti¬ 
cal works where the spin currents were found to vanish at SIN 
interfaces,^® and thus there were zero spin currents throughout 
the entire ISO coupled SINIS hybrids. Since the SINIS sys¬ 
tem considered here is in an equilibrium state, the time deriva¬ 
tive of the spin density is equal to zero.®^ Therefore, because 
the singlet superconducting electrodes considered throughout 
the paper do not support spin currents, the divergence of the 
spin current at a SIN interface is zero if the spin moment 
of the spin-active interface is zero (spin-inactive interfaces). 
This fact is clearly seen in Figs. 2(a) and 4 where the right 
SIN interface is spin-inactive. The bottom panels display the 
same components except now as a function of y along the 
junction width Jp{xo, y) at xq = 0.25^s, O.S^g, 0.75^s, and 
l.O^S. Most notably, the plots reveal a nonuniform distribu¬ 
tion of spin current densities along the junction width. From 
the plots, it is apparent that the spin current densities peak near 
the transverse vacuum boundaries of the N wire at y = 0, and 
y = Wm- Considering the top and bottom panels together, 
we conclude that the spin currents are maximally accumulated 
at the edges of the N wire near the spin-active interface and 
approximately confined within x < This in turn im¬ 

plies that the corners of the N wire near the spin-active in¬ 
terfaces possess the maximum of spin current densities.®® The 
edge accumulation of spin current densities are reminiscent 
of those previously found in nonsuperconducting mesoscale 
junctions with ISOIs.'®^^^ We emphasize that the previous 
works relied critically on externally applied magnetic and 
electric helds*®”^®’^® which can complicate the theoretical and 
experimental situations. In contrast, our hndings provide an 
alternate, simple platform which relies merely on the intrin¬ 
sic properties of the system and is devoid of any externally 


imposed conditions. Our numerical approach allows us to 
determine the precise nature of the spin and charge currents 
when varying the superconducting phase difference, (p. The 
spatial maps of charge supercurrent density (not shown) are 
constant vs position within the N wire, reflecting the charge 
conservation law. We have found that the charge supercurrent 
is governed by a sinusoidal-like current phase relation, while 
the spin currents, on the contrary, are even-functions®', i.e., 
+ 27t) = The behavior of charge supercur¬ 

rent in the low proximity limit considered here can differ from 
the ballistic regime where anomalous supercurrent-phase re¬ 
lations were found^^’®®’®'. These hndings offer an appealing 
experimental platform to examine pure edge spin currents, 
not accompanied by charge currents, solely by modulating p 
without imposing an external electromagnetic held on the sys¬ 
tem. 

A spin-active SIN interface can be ordinarily fabricated 
by coating a superconductor with a spin-active layer. 

The signatures of triplet pairings can be experimentally 
probed by means of tunneling experiments and scanning 
tunneling microscopy/spectroscopy (STM/STS) techniques^®, 
which rely on zero-energy peaks in the proximity-induced 
density of states^^’®®’®'. Technological progress allows for 
measuring high resolution spatially and energy resolved den¬ 
sity of states on a two-dimensional surface.®®’®* Therefore, the 
accumulation of triplet correlations at the boundaries or cor¬ 
ners of an N wire, and also their long-range signatures pre¬ 
dicted here, may be realized in tunneling experiments. Indeed, 
one such possibility is shown in the schematic of Fig. 1, where 
the STM tip can traverse the surface and effectively measure 
the local density of states of the entire N layer residing in the 
xy plane, producing a spatially-resolved and energy-resolved 
density of states®®’®'. Based on our hndings described thus far, 
one can expect signihcant modihcations to the local density of 
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Magnetic moment 

si={sr,sf,sn 

Sr = {0,0,0) 

s,={±sr,o,o) 

Si=i0,±Sf,0) 

si={o,o,±sn 

Rashba SOI 

jsr 

jsy 

JSZ 

jsr 

jsy 

JSZ 

jss: 

jsy 

JSZ 

Magnetic moment 

si={sr,sf,sn 

Sr = {0,0,0) 

Si={0,±sf,0) 

s,=(±sr,o,o) 

si={o,o,±sn 

Dresselhaus SOI 

jsy 

jsr 

JSZ 

jsy 

JSX 

JSZ 

jsy 

JSX 

JSZ 


TABLE I. Symmetries of the spin current components, (x, y), in 
an intrinsic spin-orbit coupled SIN IS junction with one spin-active 
interface depicted in Fig. 1. The spin moments of the left and right 
interfaces are denoted by Sw and the spin-orbit coupling is set to be 
either purely of the Rashba (a 7 ^ 0, /3 = 0) or Dresselhaus (/3 7 ^ 0, 
a = 0) type. To have succinct notation, the (x, y)-functionality of 
the spin current components is omitted in the table. The spin current 
components in similar columns have identical modulus behaviors. 

states as the STM tip moves toward the edges of the N layer 
and probes the signatures of triplet pairings, which manifest 
themselves in zero-energy peaks of the local density of states. 

Spin accumulation is a distinctive trademark of the spin- 
Hall effect.^® Therefore, the accumulation of spin current den¬ 
sities at the edges of the sample, as described in this paper, 
may be directly measurable through optical experiments such 
as Kerr rotation microscopy^®. In this scenario spatial maps 
of the spin polarizations in the entire N wire can be con¬ 
veniently imaged. A multiterminal device can also be al¬ 
ternatively employed to observe the signatures of spin cur¬ 
rents edge accumulation^When lateral leads are at¬ 
tached near the transverse vacuum edges of a two-dimensional 
S'W/5'junction (vacuum boundaries at y = 0 , and y = Wn in 
Fig. 1), the accumulated spin densities at the transverse edges 
of the N wire can result in spin current injection into the lat¬ 
eral leads, which in turn may induce a voltage drop.^'’^^’^' 

Finally, we discuss the symmetries present among the com¬ 
ponents of spin current density by varying the spin moment 
orientation of a spin-active interface, shown in Fig. 1, in a 
Rashba or Dresselhaus S/NIS hybrid. In order to systemat¬ 
ically obtain and compare results, we hrst consider a Rashba 
SOI (a 7 ^ 0, /3 = 0) and rotate Si while Sr = 0. There¬ 
after, we iterate the same procedure when the SOI is purely 
Dresselhaus (a = 0, /3 7 ^ 0). Table I summarizes the sym¬ 
metries among the components of spin current found through 
extensive numerical investigations. In the table, vector cur¬ 
rents are presented, i.e., and JP). 

The spin current components with identical spatial maps re¬ 
side in similar columns. For example, y) in the hrst 

row and column is identical to J^y{x,y) in the second row 
and hrst column. In the top row, labeled “Rashba SOI”, we 
consider Rashba SOC and rotate the spin moment of the left 
interface Si. In the bottom row (labeled “Dresselhaus SOI”), 
however, Dresselhaus SOI is considered and the spin moment 
rotations, the same as Rashba case, are iterated. As seen, when 
the moment of the spin-active interface points along the z di¬ 
rection Si = ( 0 , 0 , iSf), shows identical behaviors for 


either Rashba or Dresselhaus SOIs. However, (and 
in the presence of Rashba SOI is identical to (and J®“) 
in the presence of Dresselhaus SOI. This scenario changes 
when the moment of the spin-active interface points along the 
X or y directions. Similar symmetries to the previous case, 
i.e. Si — (0, 0,±S'f), are available, provided that we trans¬ 
form Si = (iS'f, 0,0) to Si = (0, iSf, 0), and vice a versa, 
when considering Rashba or Dresselhaus SOIs. For example, 
J®^ in the presence of Rashba SOI and Si = (±S'f, 0, 0) is 
identical to J®^ in the presence of Dresselhaus SOI, provided 
that Si = (0,±S';^,0). Under the same conditions, J®"^ (and 
J®*') in the presence of Rashba SOI is identical to J®*^ (and 
j®x) jjj jjjg presence of Dresselhaus SOI. The contents of Ta¬ 
ble I can be utilized to deduce the spin current densities in the 
presence of Rashba (Dresselhaus) SOI solely by using the re¬ 
sults of a Dresselhaus (Rashba) SOI, without any additional 
calculations. Similar transformations can take place based on 
other orientations of Si. Thus for example, this prescription 
can be used to obtain the spatial maps of spin current densities 
in a Dresselhaus S'/TV/S'junction with Si = (iS'fjOjO) and 
Sr = 0 using the data from the plots presented in Fig. 4. 


IV. CONCLUSIONS 

In conclusion, hnite-sized two-dimensional intrinsic spin- 
orbit coupled S/NIS junctions with one spin-active interface 
in the diffusive regime are theoretically studied using a quasi- 
classical approach together with spin-dependent helds obey¬ 
ing SU(2) gauge symmetries. We have computed the singlet 
and triplet correlations, and the associated spin currents in a 
S/N/S system where the interface spin moment can take ar¬ 
bitrary orientations. Using spatial maps of the singlet and 
triplet pair correlations within the two-dimensional N wire, 
we demonstrate that the combination of one spin-active inter¬ 
face and an intrinsic spin-orbit interaction (ISOI) effectively 
converts singlet pairs into long-range triplet ones. Interest¬ 
ingly, the spatial prohles illustrate that the proximity-induced 
triplet correlations are nonuniformly distributed and accumu¬ 
late at the borders of the N wire nearest the spin-active inter¬ 
face. By contrast, the spatial amplitude of the singlet correla¬ 
tions is uniform within the spin-orbit coupled N wire. The 
results suggest that the maximum singlet-triplet conversion 
takes place at the corners of the N wire nearest the spin-active 
interface. The spatial profiles of the associated spin current 
densities also demonstrate that the three components of spin 
currents accumulate the most at the edges of the N wire. Sub¬ 
sequently, the corners of the N wire near the spin-active in¬ 
terface host maximum density of spin currents.^® These re¬ 
sults are robust and independent of either the interface spin 
moment orientation or the actual type of ISOI. (We note that 
rich edge phenomena were theoretically found in hnite-size 
two-dimensional intrinsically spin orbit coupled S/FIS junc¬ 
tions in Ref. 59). We also determine the behavior of spin and 
charge currents by varying the macroscopic phase difference 
between the S banks, (f. The charge supercurrent is governed 
by the usual odd-functionality in ip, while the spin currents 
are even-functions of p, i.e. J ®((/7 + 27r) = J®(—Hence 
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by properly calibrating ip, it is possible to have pure edge spin 
currents without driving charge supercurrents. We then de¬ 
scribed experimentally relevant signatures and potential ex¬ 
periments aimed at realizing the edge phenomena predicted 
here. Our work therefore offers a simple structure consisting 
of a finite-sized intrinsic spin-orbit coupled SIN/S junction 
with one spin-active interface to generate various edge phe¬ 
nomena, such as singlet-triplet conversions, long-range prox¬ 
imity effects, and spin currents in the absence of externally 
imposed fields. 
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Appendix A: Pauli Matrices 


In Sec. II we introduced the Pauli matrices in the spin space 
and denoted them by cr = (cr^, cr^, cr^), x = 
and iz = (zz^, v^). 



We also introduced the 4x4 matrices p = {pi, p 2 , ps)' 
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Following Ref. 61, we define r'*', and po as follows; 


cr^ 0 
0 


XT' 0 

0 ai* ) ^ Po- 


CT° 0 

0 cr° ' ’ 


to unify our notation throughout the paper 7 stands for x, y, z. 
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